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A MINIATURE TWO-AXIS FLUXGATE MAGNETOMETER* 


by 

C. J. Pellerin 

Goddard Space Flight Center 
and 

M. H. Acuna 

Fairchild Hiller Corporation 


INTRODUCTION 

Ring-core magnetometers have been extensively used for magnetic measurements since their 
introduction in 1961 by W. A. Geyger (References 1 and 2). In general, the operating frequency of 
these instruments has been relatively low (5 to 10 kHz), and they have had only single-axis sensing 
capabilities because of the large number of turns required in the pickup winding. However, the 
ring-core geometry is ideally suited to multiple-axis sensing applications where the orientations 
of the sensing axes are parallel to the plane of the core. Gordon et al. (Reference 3) discussed 
the fundamental factors affecting the operation of ring-core magnetometers for different types of 
drive waveforms, and set forth criteria for selection of core materials, geometric factors, and 
operating frequency. 

This paper is an extension of that study, applied to the case of a two -axis core magnetometer 
operating at a much higher frequency and of comparatively smaller size. A linear mathematical 
model has been developed where the hysteresis loop of the core material is idealized by assuming 
constant slopes for the different regions, and is justified by comparison of the predicted and ob- 
served performance. 

This magnetometer provides two orthogonal measurements of an external magnetic field. In 
cases where the external field components are known from spherical harmonic expansions, the 
data obtained are sufficient to determine magnetic aspect as well as roll rates and approximate 
pitch angles. This is the primary application intended for the instrument. However, because of 
the high sensitivity (1 mv/y) and low construction cost, the instrument has wide application in 
geomagnetic field sensing. 

Although this paper is primarily concerned with a two -axis magnetometer, the extension to a 
three-axis magnetometer is quite simple. Laboratory prototypes have been constructed— but, for 
spinning sounding rockets, three-axis magnetometers contribute no additional information. 

^Presented at the Geoscience Electronics Symposium of the Institute of Electrical and Electronics Engineers, April 1969. 
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MAGNETOMETER CIRCUIT 


A block diagram of the instrument is shown in Figure 1. As can be observed, this is a stand- 
ard configuration for a fluxgate magnetometer. The two detectors are of the phase-sensitive type 
in order to determine the direction of the external field. In the actual instrument, the sensing core 
performs simultaneously the functions of a sensing element and of a frequency controlling device 
for the magnetic multivibrator. 

A detailed circuit diagram is shown in Figure 2. Q1 and Q2 are the switching transistors for 
the magnetic multivibrator, which drives the core periodically into saturation at a rate determined 
by the core characteristics. The second harmonic reference signal is obtained by differentiating 



Figure 1— Block diagram of two-axis ring-core magnetometer. 


Cl 



2 








the square-wave drive voltage and full-wave rectification of the resulting "spikes." T2 is an 
isolation transformer which, in addition, allows each output to be biased at any desired level, 
i.e„, 2.5 volts for telemetry applications. 

The sensitivity of each axis is adjusted by varying the damping resistors across the pickup 
windings, which are tuned to the second harmonic of the drive frequency. When measuring the 
earth’s magnetic field, this sensitivity is adjusted to 4 volts/oersted, although the intrinsic sensi- 
tivity of the instrument manufactured for sounding rocket use is of the order of 60 volts/oersted. 
This adjustment improves the linearity of the response and maintains the simplicity of the 
instrument. 


SENSING 
WINDING *1 
( 100 turns ) 


ANALYSIS OF OPERATION 

The following assumptions are made in the 
analysis: 

(a) The hysteresis loop is approximated by 
a piecewise linear function as shown in 
Figure 4. 


As the amplitude of the second harmonic reference is of the same magnitude as the supply 
voltage, the maximum dynamic range of the instrument is limited by this voltage. Typical values 
are ± 10 volts, when no bias is applied to the output. The supply voltage is regulated by means of 
an integrated -circuit voltage regulator, to 
. accommodate large variations of the input 
voltage— as are typical of sounding rocket en- 
vironments. As the operating frequency of the 
magnetic multivibrator is directly dependent 
upon its input voltage, it is important to use a 
regulator with good regulation characteristics 
and low output impedance. The frequency of 
operation of the instrument is 100 kHz, and 
intrinsic power consumption is 60 mw. 

The design of the sensing core is shown in 
Figure 3. Two orthogonal pickup windings are 
diametrically wound around the core, and occupy 
a small portion of the periphery in order to 
minimize cross -coupling effects and to improve 
angular response characteristics. The number 
of turns in each pickup winding is only about 
100; thus there is enough room to accommodate 
each winding around the core in the fashion 
shown. A plastic winding guide is used to main- 
tain the orthogonal alignment of these windings 
to within one degree. Typical core size used is 
0.5 inch in diameter and 0.125 inch in depth. 



Figure 3— Design of fhe sensing core. 
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(b) The h f e of the transistors is independent of collector current. 

(c) The drive waveform is a perfect square wave. 

(d) Stray and distributed capacitances are neglected. 

Under assumption (c), the flux variation is a linear function of time; therefore, the current 
waveform is determined by the characteristics of the hysteresis loop. As will be shown later, 
the approximation of the hysteresis loop by a piecewise linear function is justified by the fact that 
the areas of the resulting waveforms, and therefore their Fourier coefficients, are dependent only 
upon the flux- versus -time variation, which is forced to be linear by assumption (c). This is a 
good first-order approximation to the actual behavior of the circuit. The effect of a real hysteresis 
loop would be to smooth out sharp corners and spikes in the waveforms shown in Figure 4, but 
their Fourier coefficients would not be appreciably changed. 

With reference to Figures 4 and 5, we can write for the current drive function or excitation: 

0 < t < t 2 ; H(t) 

t 2 £ t < t 4 ; H(t) 

t 4 < t < T/2; H(t) = 


H + H* 

max 


1 0 a x 


(Hsat - H') (t - t 4 ) 


+ Hsat 


( 1 ) 

( 2 ) 


(H max -H sat ) 

T 72 - t 4 (t - T/2) + H m 


(3) 


where H' is given by 


H' 


H 


- 2fi , H 

' sat c 

M A 


Ufa = differential permeability, 
/LL sat = saturation permeability, 

B = flux density. 

Using the relation 


d0 dB A dH 

dt - A dt ~ dt 


4 <p 

a x 

T - 


(4) 




(5) 
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VOLTAGE AND FLUX VS TIME 


RING -CORE HYSTERESIS CURVES 
B 




Figure 5— -Hysteresis curves and waveforms. 


we can compute 


TM sat (Hmax +H ') 


4B 


( 6 ) 


where 


B m a* = M A (H sat -H c ) + ^ sat (H max -H lat ) 


(V) 


and 


Tma (H.. t -H‘) 


+ t„ 


C 4 “ 4B ' v 2 

max 


( 8 ) 
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Figure 6— The resulting waveform. 


The presence of an external magnetic field 
causes, in each half of the core, the appear- 
ance of an additional excitation H e which is 
related to the external field by means of the 
demagnetizing factor discussed later. We 
have then, from Figure 5, for the right half of 
the core, 


(h' -h +h ) 

V e max ) 


H + H' 

max 


(9) 


t 4 ( t 4 *-2) 


H . " H ' 

sat 


( 10 ) 


and analogously for the left half, 


t 


5 


T/2 + 


( H s a t +H e -H m ax)(T/2-t 4 ) 


- H 


(ID 


(»' -H.. t )(t 4 -t 2 ) 

+ Xa t - «' 


( 12 ) 


As the resulting waveform (Figure 6) is symmetrical about T/2, we can also write 


t □ t - t 

t? 1 l 5 . 


t ~ T - t 
L 8 1 t 4 


t 9 - T - t, . 


*-10 ^ t e 


t = T - t 
L n 1 *-2 • 


t = T - t 
L 12 1 * 


(13) 


The amplitudes of the resulting voltages in Figure 5 are given by 


4 NAB 


V i' 


max ( “s a t 




10 " 


(14) 
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and 


V ' = V ' 

2 ^ sat V 1 


where 

B max is given by Equation 7, and 
N = number of turns in the pickup winding, 

A = cross-sectional area of the magnetic core. 

The differential voltage output can now be decomposed into its Fourier components in order 
to determine the amplitude of the second harmonic signal as a function of the external field. The 
results of this analysis can be summarized as follows: 


f(t) - (A n cos nwt + B n sin nwt) 


where 


f(t) cos nwt dt - 0 (odd function) 


The B coefficient is then given by 


f(t) sin nwt dt 


- COS T 


XI ( 

4 77n \ 

- w ±i 

770 \ 


277nt 

COS 


3 " 4 

- - COS Y 


ltu 7 

»j» CO S 'j* 


)-4(. 


277nt 

COS J~ 


8 9 | 

- COS x I 


XL 

( 277nt 10 

2vnt u\ 

V # 
V 2 

z' 2nnt 11 

27mti 2 \ 

rm 

^COS j 

“COS J 1 

77 n 

\ cos T 

“COS 'J' J 
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For n = 2 we obtain the amplitude of the second harmonic signal, 

V 2 ( **1 ** 2 ^ , 1fn 

B 2 ~ ~2n \ cos T ~ cos T / + “* • Uoy 

Because of the complexity of the dependence of the coefficient B 2 on the external field, a closed- 
form result is meaningless; therefore, the equations described were programmed in a digital 
computer for different values and combinations of the parameters involved. The results of these 
computations are presented later on in this paper. 


THE DEMAGNETIZING FACTOR 

In order to determine the functional relationship between the external field and the excitation 
appearing in each half of the core, H e the so-called "demagnetizing factor,” must be included in 
the analysis. 

If we assume the core to be an infinitely long cylinder and solve the magnetostatic potential 
problem described by 


V 2 <D = 0 ; H = - V0> , (19) 

m m 

subject to appropriate boundary conditions, the results obtained are not satisfactory, as the pre- 
dicted excitation for a given external field is far too small and does not agree with laboratory 
measurements. 

Following the approach taken by Gordon (Reference 3), based on Bozorth’s (References 4 and 
5) work with ellipsoids and long rods, the relationship between the input signal H and the resulting 
excitation H e can be expressed as 


H e = H - H d (20) 

where H d is the self-demagnetizing field. Now 

H d = K(B-m'H) (21) 

A 

where /j.‘ is the permeability of the surrounding medium (air) and K is the demagnetizing factor. 

If 


B = mH » y-' H 
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(since y >1,000^'), then 


H d = KB = K/xH e (22) 

where y is the permeability of the core material. For long rods, the demagnetizing factor K can 
be approximated as 


K * (L/d)' 1 ’ 72 ( for y - 00 ) (23) 

where L/d is the ratio of the length of the rod to its diameter. We then can write 

. 1 . (24) 

H l + /z(L/d)~ 1 ■ 7 2 

In the case of a ring-core magnetometer, L corresponds to the length around one half of the cir- 
cumference of the core, and d to the diameter of a rod with cross-section area equivalent to that 
of the magnetic core. This assumption has been used in the analysis. If the external field is 
small, i.e. , 10" 4 oersteds, the permeability to be used in Equation 24 is the initial permeability of 
the material f x 0 . However, for larger fields this is no longer true. In this case we have approxi- 
mated the variation of the permeability with the applied field by means of an exponential function, 
i.e., 


f i 


o K 0 


(25) 


in other words, y -*■ as H p * H c ; this behavior is used to determine the value of the constant K 0 . 

As Equation 24 is now a transcendental equation, the corresponding internal field for each 
value of the assumed external field was obtained by performing an iterative calculation of Equa- 
tions 24 and 25 to the desired degree of accuracy, typically 0.1 percent. 


COMPARISON OF PREDICTED AND OBSERVED PERFORMANCE 

Figure 7 shows the predicted and actual responses of the two channels in a flight magnetometer 
with a 4-79 Mo-Permalloy core, (L/d = 46). No adjustment constants have been incorporated in 
the model to approximate actual results. The variation of the permeability with the internal field 
is assumed to be exponential, as described earlier. Material properties were obtained from 
Reference 6. 
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Figure 7— Output voltage vs external field 
(4-79 Mo-Permalloy). 



EXTERNAL MAGNETIC FIELD (oersteds) 


Figure 8— Output voltage vs magnetic field (4-79 Mo- 
Permal loy-comparison of constant and exponential yu- 0 
approximations). 


Figure 8 shows a comparison of the re- 
sults obtained assuming constant and expon- 
entially varying functions for the permeability. 
As the attenuation factor, 



1 

1 + M 0 (1/d)" 1 ' 72 


(26) 


is strongly dependent upon the permeability of 
the material, it is advantageous to use low- 
initial -permeability core materials. In the 
case of 50-50 Ni-Fe (Deltamax), /z 0 = 560; and, 
for the same L/d ratio, the attenuation factor 
is about 5 percent of the corresponding value 
for 4-79 Mo-Permalloy. In addition, the max- 
imum flux density for this material is about 
twice the value of that for 4-79 Mo-Permalloy 
—a property which, for a constant operating 
frequency, would allow the use of a core with 
one-half of the cross-section area. Hence a 
higher L/d ratio would be obtained, thus de- 
creasing even further the attenuation factor. 
However, the area of the hysteresis curve for 
50-50 Ni-Fe is about 10 times the correspond- 
ing area for 4-79 Mo-Permalloy; therefore, 
drive-power requirements are increased by 
an order of magnitude. 

The predicted response for a 50-50 Ni-Fe 
ring-core magnetometer is presented in Fig- 
ure 9, and compared to that of a 4-79 Mo- 
Permalloy instrument with same L/d ratio. 

The exponential approximation function for the 
permeability was used in both cases. 


The nature of the computer program used to analyze the equations describing the operation of 
the magnetometer allows any desired parameter to be incrementally varied, and the resulting out- 
put voltage to be displayed as a family of curves (Figure 10) for different values of the external 
field; in this case .1, .2, .3, .4, and .5 oersteds. All of these curves describe results obtained for 
4-79 Mo-Permalloy core material with represented by the exponential approximation. Curves 
for 50-50 Ni-Fe are not shown, as they are similar and experimental verification has not been 
completed. 
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EXTERNAL MAGNETIC FIELD (oersteds) 

Figure 9— Output voltage vs external field (4-79 
Mo-Permalloy and Ni-Fe core materials). 




FL - .30e 


H- = . 20e 


H e = .10e 


Figure 10— Voltage output vs input voltage 
(frequency) (4-79 Mo-Permal loy) . 


Figure 11— Voltage output vs core radius 
(50-50 Ni-Fe Permalloy). 

Figure 10 shows the output variation as a 
function of the applied input voltage. The sharp 
increase in sensitivity is due, among other 
factors, to the increase in operating frequency 
characteristic of the magnetic multivibrator. 



Figure 12— Output voltage vs differential 
permeability (4-79 Mo-Permalloy). 


The output variation as a function of the core radius is shown in Figure 11. As this radius is 
increased, the L/d ratio becomes larger, thus decreasing the attenuation factor. Figures 12, 13, 
and 14 show the effects of varying material-dependent parameters. The customary definitions are 
used, except that these parameters are defined in terms of the assumed piecewise linear B-H 
curve for the theoretical model. 
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EXTERNAL MAGNETIC FIELD ( millioersteds ) 


Figure 17 — Percentage deviation from linear fit 
vs external field (mathematical model). 

ERROR ANALYSIS 

For the flight magnetometer mentioned 
earlier, the percentage deviation of the output 
voltage from a least-squares straight-line fit 
is shown in Figures 15 and 16. The S-shaped 
curve is typical of fluxgate magnetometers, and 
is the result of the gradual decrease in sensi- 



0 50 100 150 200 250 300 350 

ANGLE (degrees) 


Figure 1 8 — Percentage deviation from cosine response 
vs angle (X-axis output, rotation in plane of core). 


tivity as the external field is increased in mag- 
nitude. The increase of as a function of the internal field and the nonlinearities of the actual 
B-H curve are probably the dominant causes of this decrease. 


Figure 17 shows a plot of the linearity error obtained from the theoretical response predicted 
by the model. As can be observed, the model satisfactorily approximates the actual error curve. 

In aspect sensing applications, it is important to determine how closely the angular response of 
the magnetometer follows a cos# type curve. Figures 18, 19, and 20 show the percentage deviation 
from an ideal cos 6 response curve as a function of the orientation angle of the instrument with 
respect to the external field. In Figures 18 and 19, the external field is rotated about an axis 
perpendicular to the plane of the core, a procedure which allows both outputs to be read out 
simultaneously; i.e., K cos 6 and K sin 6 response. In Figure 20, the external field is perpendicular 
to, and rotated about, the X- axis. 


For this analysis, the linearity error has been removed to a first-order approximation by a 
piecewise linear calibration process. The sinusoidal nature of the resulting curves is probably 
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Figure 1 9— Percentage deviation from cosine response Figure 20— Percentage deviation from cosine response 

vs angle (Y-axis output, rotation in plane of core). vs angle (X-axis output, rotation about Y-axis). 


caused by a constant angular misalignment with respect to the external field. In practice, the 
above-described errors are largely reduced during actual data reduction by performing a direct 
voltage-to-angle conversion and using the calibration curves to correct the results. 


CONCLUSIONS 

An accurate two-axis magnetometer of extreme simplicity can be constructed utilizing a 
single ring-core sensor, and sensitivities of the order of 1 mv/y can be obtained. When the ring 
core performs the dual functions of frequency controlling device and sensor element, a linear 
mathematical model can be formulated which accurately predicts observed performance. This 
model allows several parameters to be incrementally varied, and their effects on sensitivity and 
linearity of the magnetometer analyzed. 

Linearity and angular response characteristics of the actual instrument are very good and 
follow smooth, predictable patterns. This fact allows the use of correction factors during the 
data reduction process which improve the overall attainable accuracy of the instrument. 
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Appendix A 


Technical Description of the Miniature Fluxgate Magnetometer, Type FMF 


The MFM is a miniature two -axis magnetometer designed at Goddard Space Flight Center for 
rocket aspect sensing. It simultaneously measures two perpendicular components of an external 
magnetic field. Sufficient data, therefore, are obtained to determine magnetic aspect as well as 
the total field for spinning rockets. These parameters can be determined twice as often as the 
spin rate of the rocket. 

The MFM operates on the fluxgate principle. A toroidal sensor is driven into saturation at 
120 kHz. The sensor therefore radiates, but because of the toroidal geometry the resulting RFI 
is very low. When payloads are extremely sensitive to RFI, the MFM can be covered with con- 
ducting paint. This reduces the RFI to extremely low levels. 

The MFM is biased at + 2.5v and has an output which swings from 0 to 5 v in an external field 
of ± 600 millioersteds (± 60 K gammas), or approximately ± the earth’s field. Separate signal and 
power grounds are provided. The output should be loaded with 100 K ohms to obtain uniformity 
of calibration curves. 

The MFM may be secured directly to a nonmagnetic shelf with nonmagnetic hardware. Care 
must be taken to keep magnetic components away from the magnetometer (or vice versa). As the 
effect of other devices on the magnetometer falls off as l/(dist) 3 , careful placement is essential. 
Stray magnetic fields can degrade the accuracy of or saturate the magnetometer with the result 
that little or no data will be obtained, if the above precautions are not taken. Since angles are 
being measured, the sensor does not have to be on the roll axis; however, the angular relationships 
of the magnetometer with some reference point and with the solar sensor must be recorded. 

If magnetic aspect or rocket attitude is desired, a magnetometer calibration in the payload is 
an absolute necessity if accuracy is to be maintained. This requires about one-half day at the 
GSFC Magnetic Test Facility and should be one of the last things done before flight. The Test 
Facility should be reserved in advance. 

If only roll data are desired, the in-payload calibration may be omitted; but experience has 
shown that frequently a postflight decision is made to attempt to reconstruct rocket attitude with 
an uncalibrated magnetometer. This process is always difficult and sometimes impossible. 
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MFM Specifications 
(tentative) 

Input Voltage: 24v dc to 45v dc (internal over-voltage protection) 

Dimensions: 2-3/4" x li" x 1-1/8" 

Input Current: 16 milliamperes 

Range of Field: 0 ± 600 millioersteds (± 60 K gammas) 

Sensitivity: 2.5 v ± lOOmv per 600 millioersteds 
Stability of Sensitivity: ±3% 

DC Output for Zero Field: 2.5 v ± 50 mv 
Operating Temperature Range: 0°C to 60 c C 
Linearity: ±1.5% 

Frequency Response: Flat to 500 Hz 
Output Impedance: 20 K ohms 
Output Load (for standard calibration): 100 K 
Weight: 2 ounces 

Connector Pins: DEM-9-P 1) spare, 2) +28, 3) signal ground, 4) power ground, 

5) roll axis output, 6) pitch axis output, 7), 8), and 9) spares. 

Note: Signal and power ground should be externally tied together. 

FRONT TOP 
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Appendix B 


Computer Programs Used in the Analysis 


VAMC5ta]V 
V MAGS 

[I] NP + 142 

C 2 J NSEC*- 22 

[3] NS+ 100 

[4] * ENTER INPUT VOLTAGE IN VOLTS ' 

[5] EIN+U 

C 6 ] 1 ENTER EXTERNAL FIELD IN OERSTEDS ' 

C 7 ] H+- □ 

[8] 'ENTER CORE RADIUS ' 

[9] /?-□ 

[10] 0-«-3.3 

[II] HFE*- 100 

[12] P£S«-910 

[13] IMAX+EINxNSECxHFEiNPxRES 

[14] HMAX+-0.^xNP*iMAX*2xR 

[15] 'TO ENTER CORE DATA , RESPOND WITH UNDERLINED PORTION ' 

[16] 1 FOR CORE MATERIAL-SELECT 4-79 MO-PERMALLOY OR 50-50 Nl-FE ' 

[17] 2HU 

[18] -*(((Z,Z)=( *M0' , ‘A^J* ))/(]/, 0.V.0) ) 

[19] UiMODATA 

[20] ->-7 

[21] 

[22 J 7:7^4xiYPxP//m<U*£i7/xl00 0 00 0 00 

[23] (' FREQUENCY CHOSEN IN KILOHERTZ = ' ;1*7x1Q00) 

[24] LD*-( Olx((oi)*0. 5)xp) + (4x>l)*o. 5 

[25] ( ' LENGHT TO DIAMETER RATIO OF AN EQUIVALENT ROD IS ' \LD) 

[26] CRANK*-\ /( (p ,NP) Af>NSEC) ApT) ApMUDEL) ,(PNS) AoEIN) Af>Q) APR) ApHC) .(pMUSAT) ) 

[27] PHIN+MUl*-HEl+VOUT+( ( ( p , H ) , CRANK ) p 0 ) 

[28] ’70 SELECT MU-SUB-ZERO CHOOSE EITHER 1 

[29] » £0^074^7 74L0E OR POWER /1PPP0^JM/I710i7 * 

[30] Z*-((l 0 1 0)/( (Z,Z) = ( 'MO' , 'NI' ) ) ) 

[31] Zl«-E 

[32] Zl«-( (Z1,Z1)=( 'CON' , 'POW ' ) ) 

[33] Zl«-(1 0 0 1 0 0)/Zl 

[ 3 4 J Z2-*-( (+/( (Z.Zl)x( 8 4 2 1))))=(10 9 6 5 ) 

[35] 71-‘-l£ Qx(iixNSxBMAXxA*T)x( \-MUSAT + MUDF.L) 

[36] V2+VlxMUDEL*MUSAT 

[37] L+- 0 

[38] 00 : L*-L+l 

[39] -»-((Z2)/(Z3,Z4,Z5,Z6) ) 

[40] Z3 : MOCON 

[41] -+Z 7 

[42] Z 4 ; MOPOW 

[43] -*-Z 7 

[44] Z5 : NICON 

[45] ->Z 7 

[46] ZS-.NIPOW 

[47] Z7 :MUltL;l*-MU 

[48] HP+HC-(((BSAT-(MUSATx2xRC))iMUDEL)) 

[49] HEllLil+HE 

[50] 70«-4x( (MUSATx(HMAX-HSAT)) + (MUDELx(HSAT-HC) ) ) 

[51] T2*-TxMUSATx ( HMAX+HP )*TD 

[52] 74«-( TxMUDELx ( HSAT-HP) +70 ) + 72 

[53] 71<-( UiP+HMAX-HE)xT2)* (HMAX+HP) 

[54] 73+74- (HEx( 74-72 )i(HSAT-HP ) ) 

[55] TS*-(T*2) + (HSAT+HE-HMAX)x( (T*2)-TH)*(HMAX-HSAT) 

[56] 76+74+ ( HP+HE-HSAT) x ( 74-72 ) * ( HSAT-HP ) 

[57] 77+7-75 

[58] 78+7-74 

[59] 79+7-73 

[60] 710+7-76 

[61] 711+7-72 

[62] 712+7-71 

[63] ft >1// 2-*-- ( ( 71*0 2 )x( ( 10(04x74*7) )-( lo ( O 4x73 *7) ))) + (( V2 *0 2 ) x ( ( lo( 0 4x75 * 7 ) )-( lO ( 04x74 *7 ) ) ) ) 

[64] a ANZ+( ( 72*02 )x( ( lo( 0 4x2*8*7) )-( 10(0 4x2*7*7) )))+(( 71 * 02 ) x < ( 10(04x79*7) ) - ( lO( 04x78 *7 ) ) ) ) 



[6 5] n AN*i < — ( ( 71102) x( ( lo( 0 4x2*1 1*2*) )- ( lo (O 4x2*1 Or? ) ) ) ) + ( ( 72 IO 2 ) x ( ( lo (o 4 x2*l 2'l T ) ) - ( 10(04x2*1112*)) ) ) 

[66] BN l«-( ( 72102 )x( ( 20(0 4xTHT) )-( 2o(04xT2lT) ) ) ) + ( ( 71 102 ) x ( ( 20(04xT2lT) ) - ( 20 ( O 4xT6 IT ) ) ) ) 

[ 6 7 J BN2 *--( ( 71102 )x( ( 2 o(o4xT 31T) )-( 20(04xT4lT) )) ) + ( ( V2 10 2 ) x ( ( 20 ( 04xT4lT) ) -( 2o(04xT5iT) ) ) ) 

[68] BN 3«-( ( 72102 )x( ( 20(04xT7lT) )-( 20(04xT8lT) ))) + (( 71 102 ) x ( ( 2o(04xT8lT) ) - ( 20 ( O 4xT9 1 T ) ) ) ) 

[69] 5P4«~( ( 71102 )x( ( 2O(O4xT10!T) )-'( 20 ( O 4xTl 1 1 T ) ) ) ) + ( ( 72 IO 2 ) x ( ( 2o( O 4x T1 1 IT ) ) - ( 2o(04xT12lT) )) ) 

[70] •+( (T6)<(T12) )/C01 

[71] * TOTAL PERIOD EXCEEDED * 

[72] +002 

[73] CO 1 ; TP«-Ti 2 

[74] fl 4/71^( (72102 )x( (10(04 xT 21T) )-( 10(04xTllT) )))+( (71 102) x(( 1 o(04xT61T) )-(10(04xT21T) ))) 

[75] a ANT+AN 1+AN2+AN 3 +AN 4 

[76] BNT+BN1+BN2+BN3+BNH 

[77] CP<-5PT 

[7 8] ft P//Z7/[L; ] + ( _ 3o(SPTl/lPT) ) +olx ( ( ( 0 ) a( F//T*0 ) ) + ( (ANTZ 0 ) a( BNTZO ) ) ) 

[79] VOUTIL; }*-Q*CN 

[80] +iL<(pH))/CQ 

[81] C02 : ( * OUTPUT VOLTAGE * \VOUT) 

[82] * OTHER VARIABLES USED IN THIS PROGRAM ARE : 

[83J (' ELECTRICAL : »;• FFF = * ; PFF; * PFJ7F FFSJSTylflCF = * ;FFS ; * lOHMS ] •) 

[84J (• PEAK CURRENT = * \IMAX \ ' \_AMPS~\ ' } ' NUMBER OF PRIMARY TURNS = * ;PP) 

[85] (* NUMBER OF SECONDARY TURNS = * \NSEC\' Q OF PICK UP COIL = ' ;Q) 

[86] (* NUMBER OF TURNS IN PICK UP WINDING = ' \NS) 

[87] (* MAGNETIC : PEAK MAGNETIC FIELD, H = * ; HMAX ; ' [OERSTEDS] ' ) 

[88] (' PEAK INDUCTION = » ; BMAX ; ' [ GAUSS~\ * SATURATION PERMEABILITY * ; MUSAT ; * [ GA USS/OERSTED ] * ) 

[89] (' INDUCTION AT SATURATION = ' \BSAT ' \_GAUSS1 COERCIVE FORCE □ * ; HC ; » [ OERSTEDS] ' ) 

[90] (* DIFFERENTIAL PERMEABILITY = * ; MUDEL ; * l GAUSS / OERSTE Dl FIELD AT SATURATION = ' ; HSAT ; * [ OERSTEDS ] * ) 

[91] (* MU-SUB-ZERO = ’ ; MU1 [ ; 1 ] ; * [ GA USS/OERSTED'] * ) 

V 


VMOCOP[ LI ]V 

V MOCON 

[1] MU *- 14500 

[2] A TP/5’ FOP 2 PJMF// CJ/(V.P=0) FF+-F[ L ] x 2 x ( BD* 2 ) x ( 1 + ( ADiBD ) * 2 ) 1 ( MU* ( ( AD* 2 ) -BD* 2 ) ) 

[ 3 J HE+HLDH l+MUx(LD*(~l. 72) )) 

V 


V/7JC0//[[]JV 

V /VJCOP 

[1] MU*~ 5 6 0 

[2] FF+P[L]l( l+/^Ox(LO*"l. 72) ) 

V 
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[1] 

[ 2 ] 

C 3 3 

[4] 

[53 

[63 

[73 

[83 

[93 

[103 

[113 


[13 

[23 

[33 

[43 

[53 

[63 

[73 

[83 

[93 


[13 

[23 

[33 

[43 

[53 

[63 

[73 

[83 


VNIDATALU1V 
' HID AT A 

'DA%A FOR DELTAMAX AT 100 KHZ' 

BSAT+mSOQ 

HC+ 0.8 

MUDEL+ 30000 

HSAT+(BSAT+KC*MUDEL ) *A iUDEL 
MVSAT + 400 

BMAX+BSAT+ ( HMAX^HSAT ) xMUSAT 
'ENTER PEAK FLUX IN MAXWELLS' 

PHIMAX+U 

' {'CROSS SECTION AREA IS ' iPHIMAX*BMAX) 
A+PHIMAX*BMAX 


VMOPOWIQ 3V 
f MOPOW 
MU + 14500 

K+iMUDELiMU ) * 1 * 0 . 0 5 
Pl«-0 

HE+HIL 3*10 
PI 1 : HT+HE 

MU+( (PS* 0.05 )x 14500 xtf*tf£)+(ffP>0. 05 )*MUDEL 
HE+HlL]i ( l+Wi/x(iP*(“i.72) ) ) 

HE-*-{ HT+HE ) * 2 

-*•((( I (HE-HT) *HE)>0. Ol*flP)x(ffp*0) )/Pll 


7NIP0WZ 133 V 
r NIPOW 
Pl-*-0 

K+ ( MU DEL * MU ) * 1 r 0 . 3 8 
tfP«-//CL3-M0 
PI 1 : HT+HE 

MU+-( (//££0.38)x560xtf*ff£) + (tfE>0.<?8)xM£/Z?KL 
• HE+HIL1*{ 1+MU*(LD*( 1.72))) 

RE+-U1T+HE)* 2 

-►((Cl (HE-NT)*HE)>0. 01x//F)x(fl£*0) )/Pll 


V 



VMODATAIU'SV 

V MO DAT A 

ElJ 'DATA FOR 4-79 MO PERMALLOY AT 100 KHZ * 

[ 2 J BSAT + 6 0 00 

C 3 J HC+ 0.15 

[4] MUDEL + 50000 

C S j HSAT+(BSAT+HCxMUDEL ) * MUD EL 

[6J MUSAT + 100 

C 7 J BtiAX+BSAT+(HMAX-RSAT)*MVSAT 
t8j ' ENTER PEAK FLUX IK MAXWELLS' 

L 9 J PHIMAXM J 

C10J ('CROSS SECTION AREA IS * ; PHIMAXiBMAX ; ' [ CM * 2 ] * ) 
[11 3 A+PHIMAX*(BSAT*( 1 + 0.04)) 

V 


VtfyltfjSXtfCLiJV 

V SC MAG LIN X 

[1] 

[ 2 J 
C 33 

C4 J THIS IS >1 CALCULATION OF PERCENTAGE DEVIATION FROM A LEAST SQUARES FITTED STRAIGHT LINE' 
C 5 J ('SERIAL NUMBER ' ; SCI 1 ] , SCI 2 J ) 

[6] ('CHANNEL ' ;SC[3]) 

[7] 

C 8 J ' MAGNETOMETER OUTPUT VOLTAGES FOR 50£, 40* , . . . 0 . . . 40*,~50* GAMMAS RESPECTIVELY ARE' 

[9] X 

[10] Y<-( 5 43210 "1 2 “3 4 ~5)xl00 

[llj XBAR' «-+/X*p* 

[12] YBAR++/Y*pY 

[13] B++/((X-XBAR)xY)i+/((X-XBAR)x(X-XBAR)) 

[14] Sl«-1000*2 

[ 1 5 J ( 'THE NUMBER OF VOLTS PER OERSTED IS ' ;B1) 

[16] A+YBAR-B*XBAR 

[17] Z+A*B 

[18] ( ' THE ZERO OFFSET IN VOLTS IS ' ; Z ) 

[19] Z>U+Z)x£ 

[ 2 0 J P<r(Y-L)*b 

[21] AE+(+/\P)*pP 

[22] (' AVERAGE ABSOLUTE ERROR IS ‘MS; 1 PERCENT' ) 

[23] " 

[24] ' THE FOLLOWING IS A PLOT OF ERROR IN PERCENT VS THE MAGNETIC FIELD IN MI LLIOERSTEDS ' 

[25] 70 70 PLOT P VS Y 

V 
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VSKLINtQlV 

V SKLIN 

[13 ft THIS IS A PROGRAM TO FIND COS RESP ERROR TAKING LINEARITY ERR INTO ACCOUNT 

C 2 3 LIN 

C 33 Bl*-( t 24 ) - 1 

[43 Bl*-Blx 15*(O2)*360 

[53 B2+20B1 

[63 B 2 *-B 2>«500 

[73 F1+-F1 

[83 -(( |Fl[lj)>400)/Ll 

[ 9 3 B2+10B1 

[103 B 2 *-B 2x500 

[113 Ll:P+-{ (Fl)-52)*5 

C 1 2 J 'AVERAGE ABSOLUTE ERROR ' 

[133 <+/|F)*pP 

[ 1 4 J B+B ixio 

[153 70 120 PLOT P VSiB 1*360*02) 

[163 ' ANGLE IN DEGREES -+ ' 

[173 " 

[183 ( 'STANDARD VOLTAGES = ' ; X ) 

[193 ( 'STANDARD FIELD = ' \H) 

[203 {'DATA VOLTAGE = ';V) 

V 


[13 

[23 

[33 

[43 

[53 

[63 

[73 

[83 

[93 

[103 

[113 

[123 

[133 

[143 

[153 

[163 

[173 

[183 

[193 

[203 

[213 

[223 

[233 


VLJW[D3V 
LIN ; N i SPA 

ft 'THIS IS A PROGRAM TO REMOVE LINEARITY ERROR' 

Fl«-40p 0 

'ENTER THE VOLTAGES TO BE USED AS A STANDARD FOR CALIBRATION ' 


A>L3 

'ENTER THE MAGNETIC FIELDS CORRESPONDING TO THESE VOLTAGES ' 

H+U 

'ENTER THE VOLTAGES TO BE CONVERTED INTO FIELDS BASED ON 
THE ABOVE AS A STANDARD ' 

V+-U 

*«-x,*[ P X3 + UCp*3-*[( (p*)-i)3 ) 

XI*-( 2x*[13)-X[23 

X*-XI,X 

//<-tf,//[p//3 + UCp//3-//[( (ptf )-i)3 ) 

HI*-{ 2x//[l J )-//[ 23 
H+ HI , H 


N*-0 


N+-N+1 

M++/{ V[tf 3<*) 

Ml *-M + 1 

Fl[//3«-( ( VLN1-XIM] )x(;/[Ml3-//[M3 )*(*[M1]-*[M3) )+F[M3 

-*16xi(ff<plO 

'ACTUAL F MBAS' 

( 0*F1)/F1 
Fl«-< 0*F1 ) /FI 
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